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^ : 

^ i Abstract 

' We analyze the approach to chiral symmetry breaking in QCD at finite temper- 

. ature, using the functional renormalization group. We compute the running gauge 

pL| i coupling in QCD for all temperatures and scales within a simple truncated renormal- 

' izatiou flow. At finite temperature, the coupling is governed by a fixed point of the 

CN . 3-dimensional theory for scales smaller than the corresponding temperature. Chiral 

I symmetry breaking is approached if the running coupling drives the quark sector to 

^ ' criticality. We quantitatively determine the phase boundary in the plane of tempera- 

^ ■ ture and number of flavors and find good agreement with lattice results. As a generic 

^ ! and testable prediction, we observe that our underlying IR fixed-point scenario leaves 

fSJ I its imprint in the shape of the phase boundary near the critical flavor number: here, 

O ■ the scaling of the critical temperature is determined by the zero-temperature IR crit- 

Q . ical exponent of the running coupling. 

^ : 

1 Introduction and summary 

(D : 

^ • The properties of strongly interacting matter change distinctly during the transition from 

> . low to high temperatures [l], as is currently explored at heavy-ion colliders. Whereas the 
^ ! low-temperature phase can be described in terms of ordinary hadronic states, a copious 
I excitation of resonances in a hot hadronic gas eventually implies the breakdown of the 
hadronic picture; instead, a description in terms of quarks and gluons is expected to arise 
naturally owing to asymptotic freedom. In the transition region between these asymptotic 
descriptions, effective degrees of freedom, such as order parameters for the chiral or decon- 
fining phase transition, may characterize the physical properties in simple terms, i.e., with 
a simple effective action [2]. 

Recently, the notion of a strongly interacting high-temperature plasma phase has at- 
tracted much attention implying that any generic choice of degrees of freedom will not 
lead to a weakly coupled description. In fact, it is natural to expect that the low-energy 
modes of the thermal spectrum still remain strongly coupled even above the phase transi- 
tion. If so, a formulation with microscopic degrees of freedom from first principles should 
serve as the most powerful and flexible approach to a quantitative understanding of the 
system for a wide parameter range. 

In this microscopic formulation, an expansion in the coupling constant is a natural 
first step |3]. The structure of this expansion turns out to be theoretically involved [H], 
exhibiting a slow convergence behavior ^ and requiring coefficients of nonperturbative 
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origin Still, a physically well-understood computational scheme can be constructed with 
the aid of effective-field theory methods |H]. This facilitates a systematic determination 
of expansion coefficients, and the agreement with lattice simulations is often surprisingly 
good down to temperatures close to T^r [HI- The phase-transition region and the deep IR, 
however, remain inaccessible with such an expansion. 

In the present work, we use a different expansion scheme to study finite-temperature 
Yang-Mills theory and QCD in terms of microscopic variables, i.e., gluons and quarks. This 
scheme is based on a systematic and consistent operator expansion of the effective action 
which is inherently nonperturbative in the coupling. For bridging the scales from weak to 
strong coupling, we use the functional renormalization group (RG) [101 HH 1121 which is 
particularly powerful for analyzing phase transitions and critical phenomena. 

Since we do not expect that microscopic variables can answer all relevant questions in 
a simple fashion, we concentrate on two accessible problems. In the first part, we focus 
on the running of the gauge coupling driven by quantum as well as thermal fiuctuations of 
pure gluodynamics. Our findings generalize similar previous zero-temperature studies to 
arbitrary values of the temperature |13| . In the second part, we employ this result for an 
investigation of the induced quark dynamics including its back-reactions on gluodynamics, 
in order to monitor the status of chiral symmetry at finite temperature. This strategy facil- 
itates a computation of the critical temperature above which chiral symmetry is restored. 
Generalizing the system to an arbitrary number of quark fiavors, we explore the phase 
boundary in the plane of temperature and fiavor number. First results of our investigation 
have already been presented in |14| . In the present work, we detail our approach and 
generalize our findings. We also report on results for the gauge group SU(2), develop the 
formalism further for finite quark masses, and perform a stability analysis of our results. 
Moreover, we gain a simple analytical understanding of one of our most important results: 
the shape of the chiral phase boundary in the (T, A^f) plane. Whereas fermionic screening is 
the dominating mechanism for small A'f, we observe an intriguing relation between the N{ 
scaling of the critical temperature near the critical flavor number and the zero-temperature 
IR critical exponent of the running gauge coupling. This relation connects two different 
universal quantities with each other and, thus, represents a generic testable prediction of 
the phase-transition scenario, arising from our truncated RG ffow. 

In Sect. [21 we summarize the technique of RG ffow equations in the background-field 
gauge, which we use for the construction of a gauge-invariant flow. In Sect. [SI we discuss 
the details of our truncation in the gluonic sector and evaluate the running gauge coupling 
at zero and finite temperature. Quark degrees of freedom are included in Sect. [D and 
the general mechanisms of chiral quark dynamics supported by our truncated RG fiow 
is elucidated. Our findings for the chiral phase transition are presented in Sect. [HI our 
conclusions and a critical assessment of our results are given in Sect. [HI 
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2 RG flow equation in background-fleld gauge 



As an alternative to the functional-integral definition of quantum field theory, we use a 
differential formulation provided by the functional RG [101 HH 1121 • this approach, flow 
equations for general correlation functions can be constructed ^Hl- A convenient version is 
given by the flow equation for the effective average action which interpolates between 
the bare action Tk=A = S and the full quantum effective action T = Tk=o The latter 
corresponds to the generator of fully-dressed proper vertices. Aiming at gluodynamics, a 
gauge-invariant flow can be constructed with the aid of the background-fleld formalism 
[in], yielding the flow equation [T2| 

2 T^;\a A] + R,iT^;\A, A]) A 

(2) 

Here, denotes the second functional derivative with respect to the fluctuating field A, 
whereas the background-field denoted by A remains purely classical. The ghost fields are 
not displayed here and in the following for brevity, but the super-trace also includes a trace 
over the ghost sector with the corresponding minus sign. The regulator Rk in the denom- 
inator suppresses infrared (IR) modes below the scale k, and its derivative kd^Rk ensures 
ultraviolet (UV) finiteness; as a consequence, the flow of Tk is dominated by fluctuations 
with momenta ^ k'^, implementing the concept of smooth momentum-shell integrations. 

The background-fleld formalism allows for a convenient deflnition of a gauge-invariant 
effective action obtained by a gauge-flxed calculation |16| . For this, an auxiliary symmetry 
in the form of gauge-like transformations of the background held A is constructed which 
remains manifestly preserved during the calculation. Identifying the background fleld with 
the expectation value A of the fluctuating fleld at the end of the calculation, A = A, the 
quantum effective action F inherits the symmetry properties of the background fleld and 
thus is gauge invariant, T[A] = T[A, A = A]. 

The background-fleld method for flow equations has been presented in [T^: the gauge 
flxing together with the regularization lead to gauge constraints for the effective action, 
resulting in regulator-modifled Ward-Takahashi identities [IHllini, see also [213 CHI- In this 
work, we solve the flow approximately, following the strategy developed in [THlUn]- The 
property of manifest gauge invariance of the solution is still maintained by the approxi- 
mation of setting A = A already for flnite values of k. Thereby, we neglect the difference 
between the RG flows of the fluctuating and the background fleld (see [21] for a treatment 
of this difference) . The price to be paid for this approximation is that the flow is no longer 
closed [22]; i-e., information required for the next RG step is not completely provided by 
the preceding step. Moreover, this approximation satisfles some but not all constraints 
imposed by the regulator-modifled Ward-Takahashi identities (mWTI). Here we assume 
that both the information loss and the corrections due to the mWTI are quantitatively 
negligible for the flnal result. The advantage of the approximation of using Ffc[74, A = A] 
for all k is that we obtain a gauge-invariant approximate solution of the quantum theory. 
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In the present work, we optimize our truncated flow by inserting the background-field 
dependent T^"^^ into the regulator in Eq. (|T}. This adjusts the regularization to the spectral 
flow of the fluctuations 01221; it also implies a significant improvement, since larger classes 
of diagrams can be resummed in the present truncation scheme. As another advantage, the 
background-field method together with the identification A = A for all k allows to bring 
the fiow equation into a propertime form |13 | l22 t l25] which generalizes standard propertime 
fiows [211; the latter have often successfully be used for low-energy QCD models (20) • For 
this, we use a regulator Rk of the form 

X 

Rk{x)=xr{y), y:= —— , (2) 

with r{y) being a dimensionless regulator shape function of dimensionless argument. Here 
Zk denotes a wave-function renormalization. Note that both Rk and Z^ are matrix-valued 
in field space. A natural choice for the matrix entries of Z^ is given by the wave function 
renormalizations of the corresponding fields, since this establishes manifest RG invariance 
of the flow equation.^ More properties of the regulator are summarized in Appendix El 
Identifying the background field and the fiuctuation field, the fiow equation yields 

dtTk[A = A,A] = -STr9,i?,(rf )[rf = dsSTTf{s,vz)exp ( - -if) . (3) 

Here, we have introduced the (matrix-valued) anomalous dimension 

Vz ■■= -dtlnZk = -^dtZk. (4) 

The operator f{s,r]z) represents the translation of the regulator Rk into propertime space 
given by 

f{s,vz) = - riz) + {H{s) - G{s))-^dt. (5) 

The auxiliary functions on the RHS are related to the regulator shape function r{y) by 
Laplace transformation: 

h{y) = f^^, h{y) = r dsh{s)e-y\ ±H(s) = h{s), H{0) = 0, (6) 
1 + r{y) Jq ds 

3(y) = TT^y 3{y) = ds~g{s)e-y\ ±G{s) = ~g{s), G(0) = 0. (7) 

^For recent advances of an alternative approach which is based on a manifestly gauge invariant regulator, 
see [2S1- Further proposals for thermal gauge-invariant flows can be found in |2l] . 

^For the longitudinal gluon components, this implies that the matrix entry {Zk)LL is proportional to 
the inverse gauge-fixing parameter ^. As a result, this renders the truncated fiow independent of ^, and 
we can implicitly choose the Landau gauge ^ = which is known to be an RG fixed point [SI1IS21- 
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So far, we have discussed pure gauge theory. Quark fields with a mass matrix M^^ can 
similarly be treated within our framework. For this, we use a regulator of the form [2| 



i?tW = (8) 

where is a short-hand notation for ^ — ig^. Note that the quark fields live in the 
fundamental representation. This form of the fermionic regulator is chirally symmetric as 
well as invariant under background-field transformations. For later purposes, let us list the 
quark-fiuctuation contributions to the gluonic sector; the flow of Tk[A] induced by quarks 
with the regulator ^ can also be written in propertime form, 

dsTTU{s,v^,^)exp{--m'), (9) 

with [r^^"* + Rk]^^ denoting the exact (regularized) quark propagator in the background 
field. In Eq. we have introduced the anomalous dimension of the quark field, 

1]^ := -dtlnZ^. (10) 
In complete analogy to the gauge sector, we define the operator f^,{s,ri^,m) by 

Uis, T]^, rh) = ~g^{s, m)(l - r/^) + {H^{s, rh) - G^{s, rh))j-dt . (11) 

The regulator shape function r^{y) is related to the auxiliary functions appearing in the 
definition of the operator ftp{s,ri^,m) as follows 

h^{y,m) = / dsh'f'{s,m)e-y' —H^{s,m) = h'^{s,m), H^{0,m) = 0. (13) 
Jo ds 

The corresponding functions g^{y,m), g^{s,m), and G'^(s,M) are related to each other 
analogously to Eq. (jTHI). The present construction facilitates a simple inclusion of finite 
quark masses without complicating the convenient (generalized) propertime form of the 
fiow equation. 

To summarize: the functional traces in Eq. (jH} and © can now be evaluated, for 
instance, with powerful heat-kernel techniques, and all details of the regularization are 
encoded in the auxiliary functions h, g etc. Equations (jH} and ^ now serve as the starting 
point for our investigation of the gluon sector. The fiow of quark-field dependent parts of 
the effective action proceeds in a standard fashion [23, see [2H1 for reviews; in particular, 
a propertime representation is not needed for the truncation in the quark sector described 
below. 
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3 RG flow of the running coupling at finite temperature 



At first sight, the running coupling does not seem to be a useful quantity in the nonper- 
turbative domain, since it is RG-scheme and strongly definition dependent. Therefore, we 
cannot a priori associate a universal meaning to the coupling fiow, but have to use and 
interpret it always in the light of its definition and RG scheme. 

In fact, the background-field formalism provides for a simple nonperturbative definition 
of the running coupling in terms of the background-field wave function renormalization 
Zk. This is based on the nonrenormalization property of the product of coupling and 
background gauge field, gA |Tn|. The running-coupling (3g2 function is thus related to the 
anomalous dimension of the background field (cf. Eq. (fT9|) below), 

P^^,=dtg^ = {d-A + r^)g\ j] = -^dtZ^, (14) 

where we have kept the spacetime dimension d arbitrary. Since the background field can 
naturally be associated with the vacuum of gluodynamics, we may interpret our coupling 
as the response strength of the vacuum to color-charged perturbations. 



3.1 Truncated RG flow 

Owing to strong coupling, we cannot expect that low-energy gluodynamics can be described 
by a small number of gluonic operators. On the contrary, infinitely many operators become 
RG relevant and will in turn drive the running of the coupling. Following the strategy 
developed in [TH], we span a truncated space of effective action functionals by the ansatz 

r, = rriA, A] + Tf[A, A] + rf [A, A, c, c] + Tr'>, a, 4,, v]- (i5) 

Here, T^^ and T^^ represent generalized gauge-fixing and ghost contributions, which we 
assume to be well approximated by their classical form in the present work, 

rf [A, A] = 1 jjD,[A]{A-A),)\ rf [A, A, c, c] = -cD,[A]D,[A]c, D[A] = d-igA, 

(16) 

neglecting any non-trivial running in these sectors. Here, g denotes the bare coupling, 
and the gauge field lives in the adjoint representation, A^ = A^T'^, with hermitean 
gauge-group generators T'^. The gluonic part T].^ carries the desired physical infor- 
mation about the quantum theory that can be gauge-invariantly extracted in the limit 
TnA]=Tr[AA = A]. 

The quark contributions are contained in 

Tt[A,ij,ij] = I ^(i^[A] +M^^)^ + rr*[^,^], (17) 

where M^^ denotes the quark mass matrix, and the quarks transform under the funda- 
mental representation of the gauge group. The last term r^'™*[?/^, ?/;] denotes our ansatz 
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for gliion-induced quark self-interactions to be discussed in Sect. HI In Eq. (fTTj) . we have 
already set the quark wave function renormalization to = 1, which is a combined 
consequence of the Landau gauge and our later choice for r^'™*[?/^, ?/;]. 

An infinite but still tractable set of gauge-field operators is given by the nontrivial part 
of our gluonic truncation, 



(18) 



Expanding the function W(^) = Wi9 + ^W29'^ + ^W^O^ . . . , the expansion coefficients 
Wi denote an infinite set of generalized couplings. Here, Wi is identical to the desired 
background-field wave function renormalization, Zk = Wi, defining the running of the 
coupling, 

/ = k'-'z^Y, (19) 

which Eq. ()14|1 is a consequence of. This truncation corresponds to a gradient expansion 
in the field strength, neglecting higher-derivative terms and more complicated color and 
Lorentz structures. In this way, the truncation includes arbitrarily high gluonic correlators 
projected onto their small-momentum limit and onto the particular color and Lorentz 
structure arising from powers of F^. In our truncation, the running of the coupling is 
successively driven by all generalized couplings Wi. 

It is convenient to express the flow equation in terms of dimensionless renormalized 
quantities 



i9 = g^k-'^z^^e = k~Ye, 

w{^) = g^k-'m{9)^k-'Z^Ym{k'^/f). 
Inserting Eq. (fT5|l into Eq. (jS]), we obtain the flow equation for wij)): 



(20) 
(21) 



dtw = —(4 — ri)w + Mw + 



2(4vr: 



poo ( Nf 

/ -16$:5^/.n.,^)/#(.,f)/"(^W 

"^0 t i=i 5=1 



+ h{s) 



Led 



1=1 



-2f^{sw,^)f^{sw, 



N?-l 



w 



w +2i!} w' 



r]g{s) + {h{s) - g{s)) 



w 



2 E fTisw,^)f,^{swkW - f^isw,^)f^{s 
. 1=1 

2ihis) - ~g{s))^ 



w 



(w +2?9 w) 



w +2'i9 w 



' k ' 



(22) 
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where the auxiliary functions / are defined in App. El and we have used the abbreviation 
Ed = The "color magnetic" field components bi are defined by bi = |z/i|-\/2i9, where 

denotes eigenvalues of {n°-T°-) in the fundamental representation; correspondingly, bi is 
equivalently defined for the adjoint representation. Furthermore, we have used the short- 
hand notation w = w{'d) and dots denote derivatives with respect to 'd. In order to extract 
the flow equation for the running coupling, we expand the function w{'d) in powers of "(9, 



oo 

Wi 



w{d) = Y^^d\ wi = 1. (23) 

i=0 

Note that Wi is flxed to 1 by definition ((2H). Inserting this expansion into Eq. ((221), 
we obtain an infinite tower of first-order differential equations for the coefficients Wi. In 
the present work, we concentrate on the running coupling and ignore the full form of the 
function W; hence, we set — ^ for i > 2 on the RHS of the flow equation as a first 
approximation, but keep track of the fiow of all coefficients Wi. The resulting infinite tower 
of equations is of the form 

dtWi = Xi{g\ T]) + Yij{g^)dtw,, (24) 

with known functions Xj, Yij, the latter of which obeys Yij = for j > i + 1. Note that we 
have not dropped the Wi fiows, dtWi, which are a consequence of the spectral adjustment 
of the fiow. This infinite set of equations can iteratively be solved, yielding the anomalous 
dimension as an infinite power series of (for technical details, see [TTSl IHTS] ). 

2 

, = 5^a.G- with G^^^. (25) 

The coefficients be worked out analytically; they depend on the gauge group, the 

number of quark fiavors, their masses, the temperature and the regulator. Equation 
constitutes an asymptotic series, since the coefficients am grow at least factorially. This 
is no surprise, since the expansion (|23|) induces an expansion of the propertime integrals 
in Eq. (|22ll for which this is a well-understood property IHl]. A good approximation of 
the underlying finite integral representation of Eq. (f25|l can be deduced from a Borel 
resummation including only the leading asymptotic growth of the am, 



^^E"™^'"- (26) 



The leading growth coefficients are given by a sum of gluon/ghost and gluon-quark contri- 
butions. 



.iT{zd + m)r(m + 1) 



n2m 9 

ai- = 4(-2co-^ ^--^';;;7^^ ' [hL-eM)(d-2)^^r:^B2m 
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The auxiliary functions ci, C2, Zd and the moments hj, hj are defined in App. 1X1 and IbI 
The group theoretical factors and are defined and discussed in App. O The last 
term in the second line of Eq. (j?7j) contains the quark contributions to the anomalous 
dimension. The remaining terms are of gluonic origin. 

The first term in the second line has to be treated with care, since it arises from 
the Nielsen-Olesen mode in the propagator I^Hl which is unstable in the JR. This mode 
occurs in the perturbative evaluation of gradient-expanded effective actions and signals 
the instability of chromo fields with large spatial correlation. At finite temperature, this 
problem is particularly severe, since such a mode will strongly be populated by thermal 
fluctuations, typically spoiling perturbative computations |36| . 

From the flow-equation perspective, this does not cause conceptual problems, since no 
assumption on large spatial correlations of the background held is needed, in contrast to 
the perturbative gradient expansion. 

For an expansion of the flow equation about the (unknown) true vacuum state, the 
regulated propagator would be positive deflnite, + -Rfe > for /c > 0. Even without 
knowing the true vacuum state, it is therefore a viable procedure to include only the 

(2) 

positive part of the spectrum of F^ + Rk in our truncation, since it is an exact operation 
for stable background fields. At zero temperature, these considerations are redundant, 
since the unstable mode merely creates imaginary parts that can easily be separated from 
the real coupling fiow. At finite temperature, we only have to remove the unphysical 
thermal population of this mode which we do by a T-dependent regulator that screens 
the instability. As an unambiguous regularization, we include the Nielsen-Olesen mode 
for all A; > T as it is, dropping possible imaginary parts; for /c < T we remove the 
Nielsen-Olesen mode completely, thereby inhibiting its thermal excitation. Of course, a 
smeared regularization of this mode is also possible, as discussed in App. Therein, 
the regularization used here is shown to be a point of "minimum sensitivity" in a 
whole class of regulators. This supports our viewpoint that our regularization has the least 
contamination of unphysical thermal population of the Nielsen-Olesen mode. 
We outline the resummation of r] of Eq. (|26| in App. [HI yielding 

V = vt + V2+ V", (28) 

with gluonic parts ?7j^,?7^ and the quark contribution to the gluon anomalous dimension 
3 pinite integral representations of these functions are given in Eqs. (jB.lTj) . (|B.23jl . and 
(|R24jl . For pure gluodynamics, rj^ and carry the full information about the running 
coupling. 

In Fig. [U we show the result for the anomalous dimension ?7 as a function of G = ^ 
for A'c = 3 and N{ = 3 in d = 4 dimensions. For pure gluodynamics (i.e. A^f = 0), we find 
an IR stable fixed point for vanishing temperature, 

a, = [a,,8,a*,3] ~ [5.7,9.7], (29) 

^The contribution 77'' should not be confused with the quark anomalous dimension 77^ which is zero in 
our truncation. 
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Figure 1: Anomalous dimension ?7 as a function of G = ^ for Ad SU{Nc = 3) theory with 
iVf = 3 massless quark flavors at vanishing temperature. The gluonic parts rj^, rj2 and 
the quark part rj'^ contributing to the anomalous dimension rj (thick black line) are shown 
separately. The gluonic parts r]^ and agree with the results found in The figure 
shows the results from a calculation with a background field pointing into the 8-direction 
in color space. 

in agreement with the results found in [12]. The (theoretical) uncertainty is due to the 
fact that we have used a simple approximation for the exact color factors and rj^, see 
App. Ofor details. This approximation introduces an artificial dependence on the color 
direction of the background field. The extremal cases of this dependence are given by the 
3- and 8-direction in the Cartan sub-algebra, the results of which span the above interval 
for the IR fixed point. Even though this uncertainty is quantitatively large in the pure-glue 
case, it has little effect on the quantitative results for full QCD, see below. 

The inclusion of light quarks yields a lower value for the infrared fixed point a^, as can 
be seen from Fig. [H However, this lower fixed point will only be attained if quarks stay 
massless or light in the deep IR. If xSB occurs, the quarks become massive and decouple 
from the flow, such that the system is expected to approach the pure-glue fixed point. In 
any case, we can read off from Fig. ^that, already in the symmetric regime, the inclusion 
of quarks leads to a smaller coupling for scales k > k^sB, as compared to the coupling 
of a pure gluonic system. 

3.2 Running-coupling results 

For quantitative results on the running coupling, we confine ourselves to d = 4 dimensions 
and to the gauge groups SU(2) and SU(3). Of course, results for arbitrary d > 2 and other 
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Figure 2: Thermal moments as a function of for the exponential regulator. The mo- 
ments hf as well as hf are finite in the limit j ^ 0. The gluonic thermal moments hf 
grow linearly for increasing ^ due to the presence of a soft Matsubara mode, whereas the 
fermionic thermal moments hf are exponentially supressed for ^ oo. 



gauge groups can straightforwardly be obtained from our general expressions in App. iBl '^ 
To this end, a quantitative evaluation of the coupling flow requires the specification 
of the regulator shape function r(y), cf. Eq. In order to make simple contact with 
measured values of the coupling, e.g., at the Z mass or the r mass, it is advantageous to 
choose r{y) in correspondence with a regularization scheme for which the nmning of the 
coupling is sufficiently close to the standard MS running in the perturbative domain. Here, 
it is important to note that already the two-loop /?g2 coefficient depends on the regulator, 
owing to both the truncation as well as the mass-dependent regularization scheme. As 
an example, we give the two-loop l3g2 function calculated from Eq. (f22| for QCD with N^. 
colors and A^f massless quark fiavors in d = 4 dimensions: 

Pig') = - 4MN.-thtNAj^,- '--gt (30) 



(47r; 

mr^hi + NfTfht , 

-■UN^-l)Ch\-g\) + 2{H^ 

\ 2 2 



2 




45 

The moments gf^'^,h^^'^, Gf and are defined in App. El They specify the regulator 
dependence of the loop terms and depend on as is visualized in Fig.O We observe that 
even the one-loop coefficient is regulator dependent at finite temperature, but universal 
and exact at zero temperature, as it should. The latter holds, since gf^^ij = 0) = 1 



2 



and hi (-r = 0) = 1 for all admissible regulators. Using the exponential regulator, we 



^For instance, this offers a way to study nonperturbative renormalizability of QCD-like theories in extra 
dimensions as initiated in PSI for pure gauge theories. 
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find h\{^ = 0) = 2C(3), g\i^ = 0) = 1, hf'^i^ = 0) = |, = 0) = § and 

Hq{^ = 0) = C(3) for the moments at zero temperature. Using the color factors and 
T2 from App. m we compare our result to the perturbative two-loop result, 

and find good agreement to within 99% for the two-loop coefficient for SU(2) and 95% 
for SU(3) pure gauge theory. Beside this compatibility with the standard MS running the 
exponential regulator is technically and numerically convenient. 

The perturbative quality of the regulator is mandatory for a reliable estimate of abso- 
lute, i.e., dimensionful, scales of the final results. The present choice enables us to fix the 
running coupling to experimental input: as initial condition, we use the measured value of 
the coupling at the r mass scale [HHl, cts = 0.322, which by RG evolution agrees with the 
world average of as at the Z mass scale. We stress that no other parameter or scale is used 
as an input. 

The global behavior of the running coupling can be characterized in simple terms. 
Let us first concentrate on pure gluodynamics, setting iVf — > for a moment. At zero 
temperature, we rediscover the results of [131, exhibiting a standard perturbative behavior 
in the UV. In the IR, the coupling increases and approaches a stable fixed point g^ which 
is induced by a second zero of the Pg2 function, see Fig. [HI The appearance of an IR fixed 
point in Yang-Mills theories is a well-investigated phenomenon also in the Landau gauge 
[39] . Here, the IR fixed point is a consequence of a tight link between the fully dressed gluon 
and ghost propagators at low momenta which is visible in a vertex expansion jlOI. Most 
interestingly, this behavior is in accordance with the Kugo-Ojima and Gribov-Zwanziger 
confinement scenarios Even though the relation between the Landau-gauge and the 
background-gauge IR fixed point is not immediate, it is reassuring that the definition 
of the running coupling in both frameworks rests on a nonrenormalization property that 
arises from gauge invariance |42l[Tn]. Within the present mass-dependent RG scheme, the 
appearance of an IR fixed point is moreover compatible with the existence of a mass gap: 
once the scale k has dropped below the lowest physical state in the spectrum, the running 
of physically relevant couplings should freeze out, since no fluctuations are left to drive 
any further RG flow. Finally, IR fixed-point scenarios have successfully been applied also 
in phenomenological studies | ^ lii l 1 ^ 1 ^ IT7 | 1^ . 

At finite temperature, the small-coupling UV behavior remains unaffected for scales 
A; ^ T and agrees with the zero-temperature perturbative running as expected. Towards 
lower scales, the coupling increases until it develops a maximum near /c ~ T. Below, 
the coupling decreases according to a powerlaw g^ ~ k/T, see Fig. |H1 This behavior has a 
simple explanation: the wavelength of fiuctuations with momenta < is larger than the 
extent of the compactified Euclidean time direction. Hence, these modes become effectively 
3-dimensional and their limiting behavior is governed by the spatial 3d Yang-Mills theory. 
As a nontrivial result, we observe the existence of a non-Gaufiian IR fixed point also in 
the reduced 3d theory, see also Sec. 13.31 By virtue of a straightforward matching between 
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Figure 3: Running SU(3) Yang-Mills coupling aYM{k,T) as a function of k for T = 
0, 100, 500 MeV compared to the one-loop running for vanishing temperature. 



the Ad and 3d coupling, the observed powerlaw for the Ad coupling is a direct consequence 
of the strong-coupling 3d IR behavior, g'^{k ^ T) ~ g"^^^ k/T. Again, the observation of 
an IR fixed point in the 3d theory agrees with recent results in the Landau gauge [IH] . 
The 3d IR fixed point and the perturbative UV behavior already qualitatively determine 
the momentum asymptotics of the running coupling. Phenomenologically, the behavior of 
the coupling in the transition region near its maximum value is most important, which 
is quantitatively provided by the full id finite-temperature fiow equation. In addition to 
the shift of the position of the maximum with temperature, we observe a decrease of the 
maximum itself for increasing temperature. On average, the id coupling gets weaker for 
higher temperature, in agreement with naive expectations. We emphasize, however, that 
this behavior results from a nontrivial interplay of various nonperturbative contributions. 

Now, we turn to the effect of a finite number A^f of massless quark flavors. In Fig. HI we 
show the running coupling as a function of /c for T = 100 MeV and for Af = 0, . . . , 10. At 
high scales k ^ T, the running of the coupling agrees with the zero-temperature running in 
the presence of Af massless quark flavors. Towards lower scales, the coupling increases less 
strongly than the coupling of the corresponding SU(3) Yang-Mills theory, due to fermionic 
screening. At a scale ~ T, the coupling reaches its maximum. Below this scale, the 
quarks decouple from the flow, since they only have hard Matsubara modes and, hence, 
the coupling universally approaches the result for pure Yang-Mills theory. Furthermore, 
we observe that, for an increasing number of quark flavors, the maximum of the coupling 
becomes smaller and moves towards lower scales. Both effects are due to the fact that the 
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Figure 4: Running SU(3) coupling as{k, T) as a function of /c for T = 100 MeV for different 
number of quark flavors A'f = 0, 1, 2, . . . , 10 (from top to bottom). For k <^T, the coupling 
shows universal behavior, owing to the attraction of the pure- glue IR flxed point. 

anomalous dimension rj becomes smaller for an increasing number of quark flavors. 

Again, we stress that the results for the coupling with dynamical quarks have not yet 
accounted for xSB, where the quarks become massive and decouple from the flow. This 
will be discussed in the following sections. For temperatures or flavor numbers larger than 
the corresponding critical value for xSB, our results so far should be trustworthy on all 
scales. 

3.3 Dimensionally reduced high-temperature limit 

As discussed above, the running coupling for scales much lower than the temperature, 
k <^T, is governed by the IR flxed point of the 3-dimensional theory. More quantitatively, 
we observe that the flow of the coupling is completely determined by r/j^ for |^ ^ 1; the 
quark contributions decouple from the flow in this limit, since they do not have a soft 
Matsubara mode. Therefore, we flnd an IR flxed point at flnite temperature for the id 
theory at = 0. In the limit |^ ^ 1, the anomalous dimension Eq. (f28|) is given by 



V{T » k) ^ <(T » k) =: vTig', B = 73d {^9')\ 



(32) 



where 73^ is a number which depends on A',,: 



32C(|)(i-2V2)r(|)r(| + ^r)^ 
(47r)4r(|)r(zr + 1) 




(33) 
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We refer to App. [Blfor the definition of the constants zf" and cf. In the high-temperature 
hmit, we can solve the differential equation (fT^ for analytically, 

9\ =■■ aim = --W^ -TT ^ 737 ^ + Oim. (34) 

-»i (73rf(i)4 - const.)^ ^ 

The RHS explains the shape of the running coupling for small k/T in Fig. El The factor 

_4 

73/ is the fixed point value of the dimensionless 3d coupling (73^, as can be seen from its 
relation to the dimensionless coupling g"^ in four dimensions: 

2 -^2 2^2 /oc\ 

93d -=-^9 ^9 = 7^ 93d- (35) 

Comparing the right-hand side of Eq. (f3i|) and (f35|) , we find that the fixed point for Nf. = ?) 
in three dimensions is given by: 

^ ^ = [al%, al%] ^ [2.70, 2.77] (36) 

Again, the uncertainty arises from our ignorance of the exact color factors r^, see App. |B1 
and App. O 

On the other hand, the fixed point of the 3d theory is determined by the zero of the 
corresponding (3 function. In fact, rff^g"^^ |^) is identical to the 3d anomalous dimension 
Vsd{g3d)j as can be deduced from the pure 3d theory, and we obtain 

dtih') = 9tgl, = {m{9ld) - 1)^?L, (37) 

as suggested by Eq. (|THl . Since 773^ is a monotonously increasing function, we find a 3d IR 
fixed point for g1^ ^ = 73/ which coincides with the result above. 



4 Chiral quark dynamics 

Dynamical quarks infiuence the QCD fiow by two qualitatively different mechanisms. First, 
quark fiuctuations directly modify the running coupling as already discussed above; the 
nonperturbative contribution in the form of i]'^ in Eq. accounts for the screening nature 
of fermionic fiuctuations, generalizing the tendency that is already visible in perturbation 
theory. Second, gluon exchange between quarks induces quark self-interactions which can 
become relevant in the strongly coupled IR. Both the quark and the gluon sector feed 
back onto each other in an involved nonlinear fashion. In general, these nonlinearities have 
to be taken into account and are provided by the fiow equation. However, we will argue 
that some intricate nonlinearities drop out or are negligible for locating the chiral phase 
boundary in a first approximation. 
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Working solely in d = 4 from here on, let us now specify the last part of our truncation: 
the effective action of quark self-interactions r^'™*[?/;, ?/;], introduced in Eq. (fT7|l . In a 
consistent and systematic operator expansion, the lowest nontrivial order is given by [52] 



A_(V-A) + A+(V+A) + A,(S-P) + AvA[2(V-A)^^j+ (l/Ar,)(V-A)] . (38) 



The four-fermion interactions occurring here have been classified according to their color 
and flavor structure. Color and flavor singlets are 

(V-A) = (^A7m^)' + (V'TmTs^)', (39) 
(V+A) = (^^^^)2 _ (^^^^5^)2, (40) 

where (fundamental) color . . . ) and flavor (x, ^, . . . ) indices are contracted pairwise, 
e.g., (ipip) = (^f^f). The remaining operators have non-singlet color or flavor structure, 

(s-p) = ir^^f - iri^^^? = i^M? - (#75^^.^)', 

(V-A)^^j = (V'7m^»2 + (^7^75T'^^/;)2, ^4^^ 

where {tp^tl^^y = ijj^ip^ijj^ip^, etc., and {T"-)ij denotes the generators of the gauge group 
in the fundamental representation. The set of fermionic self-interactions introduced in 
Eq. (I38jl forms a complete basis. Any other pointlike four-fermion interaction which is 
invariant under SU(A'^(;) gauge symmetry and SU(A^f)L x SU(A^f)R flavor symmetry is re- 
ducible by means of Fierz transformations. UA(l)-violating interactions are neglected, 
since we expect them to become relevant only inside the xSB regime or for small iVf; since 
the lowest-order UA(l)-violating term schematically is ~ {ipip)^', larger iVf correspond to 
larger RG "irrelevance" by naive power-counting. For A^f = 1, such a term is, of course, im- 
portant, since it provides for a direct fermion mass term; in this case, the chiral transition 
is expected to be a crossover. Dropping the UA(l)-violating interactions, we thus conflne 
ourselves to Af > 2. 

We emphasize that the A's are not considered as independent external parameters as, 
e.g., in the Nambu-Jona-Lasinio model. More precisely, we impose the boundary condition 
Aj ^ for A; ^ A — i> 00 which guarantees that the A's at k < A are solely generated by 
quark-gluon dynamics, e.g., by IPI "box" diagrams with 2-gluon exchange. 

As a severe approximation, we drop any nontrivial momentum dependencies of the 
A's and study these couplings in the point-like limit X{\pi\ <^ k). This inhibits a study 
of QCD properties in the chirally broken regime, since mesons, for instance, manifest 
themselves as momentum singularities in the A's. Nevertheless, the point-like truncation 
can be a reasonable approximation in the chirally symmetric regime; this has recently been 
quantitatively conflrmed for the zero-temperature chiral phase transition in many-flavor 
QCD [nni, where the regulator independence of universal quantities has been shown to 
hold remarkably well even in this restrictive truncation. By adopting the same system at 
finite T, we base our truncation on the assumption that quark dynamics both near the 
finite-T phase boundary as well as near the many-fiavor phase boundary [HI] is driven by 
qualitatively similar mechanisms. 
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The resulting flow equations for the A's are a straightforward generalization of those 
derived and analyzed in [HSl 1^0] to the case of finite temperature. Introducing the dimen- 
sionless renormalized couplings 

A. = k^Xi, (42) 



(recall that Z^, = 1 in our truncation), the flows of the quark interactions read 



dfX^ 



dtX^ 



2A_ - AvJ\'l 



(FB) 



■g'x^ - Sg'X 



VA 



^ /(FB) 



12 + 9N^ 



(43) 



-8v/^^\^ - NiN.iXl + A^^) + A^ - 2(iV, + iVf)A_AvA + iVfA+A, + 2A^a}, 



(FB) 



2X+-AvJ\'; 



-—g^x, 

Nr 



-v4\ 



(FB) 



12 + 3iV2 ■ 



(44) 



-Sv/P\^ - 3A^ - 2N^NiX.X+ - 2A+(A_ + {N^ + A^f)AvA) + iVfA.A^ 
+AvAAcr + jAo-^j, 



dtX„ = 2X„- Av4ir [6C2(iVc) g^X, - Qg^X+] - -v^lf^, 
-8t;4/r)|2iV,A2-2A_A,-2iVfA,AvA-6A+A,}, 



24 - 9Ar2 



' g' 



OA 



t^VA 



2AvA-4t;44?^ 



■fif^AvA - 35'^A- 



,(FB) 



24-3A^2 ■ 

-g 



(45) 



(46) 



-8v/P{ - (TVc + iVf ) + 4A- A VA - iN.Xl } . 



Here, C2{Nc) = {N^ — l)/{2Nc) is a Casimir operator of the gauge group, and V4 = l/(327r^). 
For better readability, we have written all gauge-coupling-dependent terms in square brack- 
ets, whereas fermionic self-interactions are grouped inside braces. The threshold functions 
l[ 1 n,i depend on the details of the regularization, see App. El for zero quark mass 

and vanishing temperature, these functions reduce to simple positive numbers, see, e.g., 
Eqs. (IA.26|1 and (jA.29|l .^ For quark masses and temperature becoming larger than the 



regulator scale /c, these functions approach zero, which reflects the decoupling of massive 
modes from the flow. 

Within this set of degrees of freedom, a simple picture for the chiral dynamics arises: for 
vanishing gauge coupling, the flow is solved by vanishing Aj's, which defines the Gauftian 
fixed point. This fixed point is IR attractive, implying that these self-interactions are RG 
irrelevant for sufficiently small bare couplings, as they should. At weak gauge coupling, the 
RG fiow generates quark self-interactions of order A ~ (?^, as expected for a perturbative 
IPI scattering amplitude. The back-reaction of these self-interactions on the total RG 
fiow is negligible at weak coupling. If the gauge coupling in the IR remains smaller than 



^Here, we ignore a weak dependence of the threshold functions on the anomalous quark and gluon 
dimensions which were shown to influence the quantitative results for the present system only on the 



percent level, if at all 
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Figure 5: Sketch of a typical (3 function for the fermionic self-interactions Aj: at zero gauge 
coupling, (7 = (upper solid curve), the Gaufiiian fixed point Aj = is IR attractive. For 
small > (middle/blue solid curve), the fixed-point positions are shifted on the order 
of (7^. For gauge couplings larger than the critical coupling g > (lower/green solid 
curve), no fixed points remain and the self-interactions quickly grow large, signaling xSB. 
For increasing temperature, the parabolas become broader and higher, owing to thermal 
fermion masses; this is indicated by the dashed/red line. 

a critical value g < gcr, the self-interactions remain bounded, approaching fixed points 
in the IR. These fixed points can simply be viewed as order-r?^ shifted versions of the 
Gaufiiian fixed point, being modified by the gauge dynamics. At these fixed points, the 
fermionic subsystem remains in the chirally invariant phase which is indeed realized at high 
temperature. 

If the gauge coupling increases beyond the critical coupling g > gen the above-mentioned 
IR fixed points are destabilized and the quark self-interactions become critical. This can 
be visualized by the fact that dtXi as a function of Aj is an everted parabola, see Fig. El 
for g = gcT, the parabola is pushed below the Aj axis, such that the (shifted) Gauftian fixed 
point annihilates with the second zero of the parabola. In this case, the gauge-fiuctuation- 
induced A's have become strong enough to contribute as relevant operators to the RG flow. 
These couplings now increase rapidly, approaching a divergence at a finite scale k = k^sB- 
In fact, this seeming Landau-pole behavior indicates xSB and, more specifically, the for- 
mation of chiral condensates. This is because the A's are proportional to the inverse mass 
parameter of a Ginzburg-Landau effective potential for the order parameter in a (partially) 
bosonized formulation, A ~ Thus, the scale at which the self-interactions formally 

diverge in our truncation is a good measure for the scale k^sB where the effective potential 
for the chiral order parameter becomes flat and is about to develop a nonzero vacuum 
expectation value. 

Whether or not chiral symmetry is preserved by the ground state therefore depends on 
the coupling strength of the system, more specifically, the value of the gauge coupling g 
relative to the critical coupling g^r which is required to trigger xSB. Incidentally, the critical 
coupling gcr itself can be determined by algebraically solving the fixed-point equations 
dt\{K) = for that value of the coupling, g = gcr, where the shifted Gaufi.ian fixed point 
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is annihilated. For instance, at zero temperature, the SU(3) critical coupling for the quarks 
system is = glJ{^T^) — 0.8 being only weakly dependent on the number of flavors 
[50].^ Since the IR fixed point for the gauge coupling is much larger a^, > ctcr (for not too 
many massless flavors), the QCD vacuum is characterized by xSB. The same qualitative 
observations have already been made in [HI] in a similar though smaller truncation. The 
existence of such a critical coupling also is a well-studied phenomenon in Dyson-Schwinger 
equations |55| . 

As soon as the the quark sector approaches criticality, also its back-reaction onto the 
gluon sector becomes sizable. Here, a subtlety of the present formalism becomes impor- 
tant: identifying the fluctuation fleld with the background held under the flow, our approx- 
imation generally does not distinguish between the flow of the background-field coupling 
and that of the fluctuation-field coupling. In our truncation, differences arise from the 
quark self-interactions. Whereas the running of the background-fleld coupling is always 
given by Eq. , the quark self-interactions can contribute directly to the running of the 
fluctuation-field coupling in the form of a "vertex correction" to the quark-gluon vertex. 
Since the fluctuation-field coupling is responsible for inducing quark self-interactions, this 
difference may become important. In |52| . the relevant terms have been derived with the 
aid of a regulator-dependent Ward-Takahashi identity. The result hence implements an 
important gauge constraint, leading us to 

dtg' = V9'-4^v/p 9,VcA, (47) 

1 - 2vJ\ ' }2 Ci\ 
= 1 + A^f, c+ = 0, c_ = -2, cvA = -2iVf, 

with rj provided by Eq. (|28|1 in our approximation. In principle, the approach to xSB can 
now be studied by solving the coupled system of Eqs. (|47|) . and ()43|l - (|46|) . However, a 
simpler and, for our purposes, sufficient estimate is provided by the following argument: 
if the system ends up in the chirally symmetric phase, the Aj's always stay close to the 
shifted Gaufiian fixed point discussed above; apart from a slight variation of this fixed- 
point position with increasing g^, the dtXi flow is small and vanishes in the IR, dtXi 0. 
Therefore, the additional terms in Eq. (jlTIl are negligible for all k and drop out in the 
IR. As a result, the behavior of the running coupling in the chirally symmetric phase is 
basically determined by t] alone, as discussed in the preceding section. In other words, 
the difference between the fluctuation-field coupling and the background-field coupling 
automatically switches off in the deep IR in the symmetric phase in our truncation. 

Therefore, if the coupling as predicted by Pg2 ^ rjg'^ alone never increases beyond the 
critical value g^^ for any k, the system is in the chirally symmetric phase. In this case, it 
suffices to solve the g'^ flow and compare it with g'^^ which can be deduced from a purely 
algebraic solution of the fixed-point equations, dtXi{X^) = 0. 

^Of course, the critical coupling is a non-universal value depending on the regularization scheme; the 
value given here for illustration holds for a class of regulators in the functional RG scheme that includes 
the most widely used linear ("optimized") and exponential regulators. 
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If the coupling as predicted by (3g2 ~ rig"^ alone approaches for some finite scale /Ccr, 
the quark sector becomes critical and all couplings start to flow rapidly. To the present 
level of accuracy, this serves as an indication for xSB. Of course, if the gauge coupling 
dropped quickly for decreasing A;, the quark sector could, in principle, become subcritical 
again. However, this might happen only for a marginal range of ~ if at all. For even 
larger gauge coupling, the flow towards xSB is unavoidable. 

Inside the xSB regime, also the induced quark masses back-react onto the gluonic flow 
in the form of a decoupling of the quark fluctuations, i.e., rj^ in Eq. (f28|) approaches zero. 
However, the present truncation does not allow to explore the properties of the xSB sector; 
for this, the introduction of effective mesonic degrees of freedom along the lines of [SBl IHE] 
is most useful and will be employed in future work. 

5 Chiral phase transition 

Let us now discuss our results for the chiral phase transition in the framework presented 
so far. As elucidated in the previous section, the breaking of chiral-symmetry is triggered 
if the gauge coupling g"^ increases beyond g^^, signaling criticality of the quark sector. We 
study the dependence of the chiral symmetry status on two parameters: temperature T 
and number of (massless) flavors iVf. As already discussed in Sect. El the increase of the 
running coupling in the IR is weakened on average for both larger T and larger A^f. In 
addition, also g^ depends on T and Af, even though the Af dependence is rather weak. 

The T dependence of g^^ has a physical interpretation: at flnite T, all quark modes 
acquire thermal masses which leads to a quark decoupling for k < T. Hence, stronger in- 
teractions are required to excite critical quark dynamics. Technically, this T/k dependence 
is a direct consequence of the T/k dependence of the threshold functions /j^^2^\ /[™^ in 
Eqs. (HTH) - see App. El for their definition. Since the threshold functions decrease 
with increasing temperature, the Aj parabolas visualized in Fig. become broader with 
a higher maximum; hence, the annihilation of the Gauftian fixed point by pushing the 
parabola below the Aj axis requires a larger 

At zero temperature and for small Af, the IR fixed point of the running coupling is 
far larger than g"^^, hence the QCD vacuum is in the xSB phase. For increasing T, the 
temperature dependence of the coupling and that of g"^^ compete with each other. This is 
illustrated in Fig. El where we show the running coupling = |^ and its critical value 

2 

ctcr = ^ for T = 130 MeV and T = 220 MeV as a function of the regulator scale k. The 
intersection point /Ccr between both marks the scale where the quark dynamics becomes 
critical. Below the scale k^r, the system runs quickly into the xSB regime. We estimate the 
critical temperature T^r as the lowest temperature for which no intersection point between 
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Figure 6: Running QCD coupling as{k,T) for A^f = 3 massless quark flavors and A'c = 
3 colors and the critical value of the running coupling acr{k,T) as a function of k for 
T = 130 MeV (left panel) and T = 220 MeV (right panel). The existence of the (as,acr) 
intersection point in the left panel indicates that the xSB quark dynamics can become 
critical for T = 130 MeV. 



as and Ocr occurs.^ We find 

r„ ^ 186 ± 33 MeV for A^f = 2, 

T„ ^ 161 ±31 MeV for A^f = 3, (48) 

for massless quark flavors in good agreement with lattice simulations The errors 

arise from the experimental uncertainties on ag |38| . The theoretical error owing to the 
color-factor uncertainty turns out to be subdominant by far, see Fig. El Dimensionless 
observable ratios are less contaminated by this uncertainty of ttg. For instance, the relative 
difference for T^r for Nf = 2 and 3 flavors is 



T,^f=2 - T^f=3 +0.018 
" " - 0.144 , (49) 



(T,f=2 + Tc^f=3)/2 ■ -0.013' 



in reasonable agreement with the lattice value of ~ 0.12 ^ 

For the case of many massless quark flavors A'f, the critical temperature is plotted in 
Fig. m We observe an almost linear decrease of the critical temperature for increasing A^f 
with a slope of AT^r = T{Nf) —T{Nf+l) ^ 25 MeV. In addition, we find a critical number 
of quark flavors, N^^ ^ 12.9, above which no chiral phase transition occurs. This result 

^Strictly speaking, this simplified analysis yields a sufficient but not a necessary criterion for chiral- 
symmetry restoration. In this sense, our estimate for is an upper bound for the true Tcr- Small 
corrections to this estimate could arise, if the quark dynamics becomes uncritical again by a strong decrease 
of the gauge coupling towards the IR, as discussed in the preceding section. 

^Even this comparison is potentially contaminated by fixing the two theories with different fiavor content 
in different ways. Whereas lattice simulations generically keep the string tension fixed, we determine all 
scales by fixing a at the r mass scale, cf. the discussion below. 
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Figure 7: Chiral-phase-transition temperature Tcr versus the number of massless quark 
flavors A'f for A'f > 2. The flattening at A'f > 10 is a consequence of the IR flxed-point 
structure. The dotted line depicts the analytic estimate near A^^" which follows from the 
flxed-point scenario (cf. Eq. (f55|) below). Squares and triangles correspond to calculations 
with a background held in the 8- and 3-direction of the Cartan, respectively. The theoretical 
uncertainty which is given by the difference between both is obviously negligible in full 
QCD. 

for N^^ agrees with other studies based on the 2-loop P function [SI]. However, the precise 
value of N^^ has to be taken with care: for instance, in a perturbative framework, N^^ 
is sensitive to the 3-loop coefflcient which can bring A^" down to A^^'' ~ 10 [S]]. In our 
nonperturbative approach, the truncation error can induce similar uncertainties; in fact, 
it is reassuring that our prediction for A^^ lies in the same ball park as the perturbative 
estimates, even though the details of the corresponding Pg2 are very different. This suggests 
that our truncation error for A^^'' is also of order 0(1). We expect that a more reliable 
estimate can be obtained even within our truncation by a regulator [Hl^ fTK]. 

A remarkable feature of the T, A^f phase diagram of Fig. U\ is the shape of the phase 
boundary, in particular, the flattening near A^". In fact, this shape can be understood 
analytically, revealing a direct connection between two universal quantities: the phase 
boundary and the IR critical exponent of the running coupling. 

Before we outline the argument in detail, let us start with an important caveat: vary- 
ing A^f unlike varying T corresponds to an unphysical deformation of a physical system. 
Whereas the deformation itself is, of course, unambiguously deflned, the comparison of 
the physical theory with the deformed theory (or between two deformed theories) is not 
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unique. A meaningful comparison requires to identify one parameter or one scale in both 
theories. In our case, we always keep the running coupling at the r mass scale fixed to 
a(mT-) = 0.322. Obviously, the couplings in the two theories are different on all other 
scales, as are generally all dimensionful quantities such as Aqcd- There is, of course, no 
generic choice for fixing the corresponding theories relative to each other. Nevertheless, we 
believe that our choice is particularly useful, since the r mass scale is close to the transition 
between perturbative and nonperturbative regimes. In this sense, a meaningful comparison 
between the theories can be made in both regimes, without being too much afflicted by 
the choice of the fixing condition. 

Let us now study the shape of the phase boundary for small N^. Once the coupling is 
fixed to a{mr) = 0.322, no free parameter is left. As a crude approximation, the mass scale 
of all dimensionful IR observables such as the critical temperature T^r is set by the scale k^o 
where the running gauge coupling undergoes the crossover from small to nonperturbatively 
large couplings (for instance, one can define the crossover scale k^o from the inflection point 
of the running coupling in Fig. As an even cruder estimate, let us approximate kco by 
the position of the Landau pole of the perturbative one-loop running coupling.^ The latter 
can be derived from the one-loop relation 

^ = ^^ + 47r6oln— , bo = — ( -N, - -NA . 50 

a[k) a[mr) rrir o-n \ 6 ^ J 

Deflning /Cco by the Landau-pole scale, 1/q;(A;co) = 0, and estimating the order of the critical 
temperature by T„ ~ /cco, we obtain 

~ Q-^^b^l^var) e" "'^c'^Cmr) (l _ giVf + C)((eA^f)^)) , (51) 

where e = YUN^^aim ) — 0.107 for = 3. This simple estimate hence predicts a linear 
decrease of the phase boundary Tcr{N{) for small A'f, as is conflrmed by the full solution 
plotted in Fig. [3 Actually, this estimate is also quantitatively accurate, since it predicts 
a relative difference for for Ni=2 and 3 flavors of A ~ 0.146 which is in very good 
agreement with the full result, given in Eq. (|i9|) . We conclude that the shape of the phase 
boundary for small A^f is basically dominated by fermionic screening. 

For larger Af, the above estimate can no longer be used, because neither one-loop 
perturbation theory nor the Af expansion are justifled. For values of Af close to the 
critical value A^", a different analytic argument can be made: here the running coupling 
has to come close to its maximal value in order to be strong enough to trigger xSB. The 
maximal value is, of course, close to the IR flxed point value attained for T = 0. Even 
though at flnite T the coupling is eventually governed by the ?>d flxed point implying a 
linear decrease with k, the xSB properties will still be dictated by the maximum coupling 
value, which roughly corresponds to the T = flxed point. In the flxed-point regime, we 
can approximate the (3g2 function by a linear expansion about the flxed-point value, 

^ dtg' = -e (/ - gl) + 0((/ - glf), (52) 

^Actually, this is a reasonable estimate, since the N( dependence of fcco, which is all that matters in the 
following, is close to the perturbative behavior. 
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6.39 


5.50 


4.99 


4.41 


3.82 


3.19 


2.58 


1.97 


1.42 


0.95 


0.57 



Table 1: The "critical exponent" G for different values of A^f for T = 0. 

where the universal "critical exponent" 9 denotes the (negative) first expansion coefficient. 
We know that 9 < 0, since the fixed point is IR attractive. For vanishing temperature, we 
find an approximate linear dependence of 9 on A'f, cf. Tab. QJ 

The solution of Eq. (l52|l for the running coupling in the fixed-point regime reads 

^'•-(^)'' '''' 

where the scale /cq is implicitly defined by a suitable initial condition (to be set in the fixed- 
point regime) and is kept fixed in the following. It provides for all dimensionful scales in 
the following and is related to the initial r mass scale by RG evolution. Our criterion for 
xSB to occur is that g'^{k) should exceed g"^^ for some value of /c = kcr- We expect that 
this scale kcr is generically somewhat larger than the temperature, since for k smaller than 
T the coupling decreases again owing to the 3d fixed point. This allows us to ignore 
the T dependence of the running coupling and of the critical coupling as a rough 
approximation, since the T dependence of the threshold functions is rather weak for T <k. 
From Eq. (f53|) and the condition g^{kcr) = fi'cn derive the estimate 

kcr^ko{gl-gly^. (54) 

This scale k^r plays the same role as the crossover scale kco in the small- A^f argument given 
above: it sets the scale for T^r ~ kcr, with a proportionality coefficient provided by the 
solution of the full fiow. To conclude the argument, we note that the IR fixed-point value 
gl roughly depends linearly on A^f, since the quark contribution to the coupling fiow r]'^ is 
linear in A^f. From Eq. (f5i|) . we thus find the relation 

Tc.~A;o|A^f-Arf"|-i, (55) 

which is expected to hold near A^^"" for A'f < N^"^. Here, 9 should be evaluated at A^".^^ 
Relation (|55|l is an analytic prediction for the shape of the chiral phase boundary in the 
(T, A'f) plane of QCD. Remarkably, it relates two universal quantities with each other: the 
phase boundary and the IR critical exponent. 

This relation can be checked with a fit of the full numerical result parametrized by the 
RHS of Eq. |55|l . In fact, the fit result, 9fit — —0.60, determined from the phase boundary 

^"Indeed, this assumption is justified, since we find in the full calculation that ka » T for large iVf and 
for temperatures in the vicinity of the critical temperature . 

Accounting for the N( dependence of Q by an expansion around N^ '^ yields mild logarithmic corrections 
to Eq. (EH). 
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agrees with the direct determination of the critical exponent from the zero-temperature (3 
function, Q{N^^ ~ 12.9) ~ —0.60, within a one-percent accuracy (cf. Table H)). The fit is 
depicted by the dashed line in Fig. [3 In particular, the fact that |9| < 1 near A^" explains 
the flattening of the phase boundary near the critical flavor number. 

Qualitatively, relation is a consequence of the IR fixed-point scenario predicted 
by our truncated flow equation. We emphasize, however, that the quantitative results for 
universal quantities such as 9 are likely to be affected by truncation errors. These can 
be reduced by an optimization of the present flow; we expect from preliminary regulator 
studies that more reliable estimates of yield smaller absolute values and, thus, a more 
pronounced flattening of the phase boundary. 

We are aware of the fact that the relation (|55|) is difficult to test, for instance, by lattice 
gauge theory: neither the fixed-point scenario in the deep IR nor large flavor numbers 
are easily accessible, even though there are promising investigations that have collected 
evidence for the IR fixed-point scenario in the Landau gauge [SHI ED] (see also (HU IS21 EBl) 
as well as the existence of a critical fiavor number [HI]. Given the conceptual simplicity 
of the fixed-point scenario in combination with xSB, further lattice studies are certainly 
worthwhile. 

6 Conclusions and outlook 

We have obtained new nonperturbative results for the chiral phase boundary of QCD in 
the plane spanned by temperature and quark fiavor number. Our work is based on the 
functional RG which provides for a functional differential formulation of QCD in terms of 
a flow equation for the effective action. We have studied this effective action from flrst 
principles in a systematic and consistent operator expansion which is partly reminiscent to 
a gradient expansion. We consider the truncated expansion as a minimal approximation of 
the eff'ective action that is capable to access the nonperturbative IR domain and address 
the phenomenon of chiral symmetry breaking. 

In the gluon sector, this truncation provides for a stable flow of the gauge coupling, 
running into a fixed point in the IR at zero temperature in agreement with the results of 
[in] for the pure glue sector. As a new result, we find that the Sc? analogue of this IR 
fixed point governs the fiow of the gauge coupling at finite temperature for scales k <^T. 
Our truncation in the quark sector facilitates a description of critical dynamics with a 
gluon-driven approach to xSB. The resulting picture for xSB is comparatively simple: 
xSB requires the coupling to exceed a critical value ga- Whether or not this critical 
coupling is reached depends on the RG fiow of the gauge coupling. The IR fixed-point 
scenario generically puts an upper bound on the maximal coupling value which depends 
on the external parameters such as temperature and quark fiavor number. Of course, the 
interplay between the gluon and quark sectors in general, and between gauge coupling and 
critical coupling in particular, is highly nonlinear, since both sectors back-react onto each 
other in a manner which is quantitatively captured by the fiow equation. 

The resulting phase boundary in the (T, iVf) plane exhibits a characteristic shape which 
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can analytically be understood in terms of simple physical mechanisms: for small A^f, 
we observe a linear decrease of as a function of A'f; this is a direct consequence of 
the charge-screening properties of light fermions. Also, this screening nature is ultimately 
responsible for the existence of a critical flavor number iV" above which the system remains 
in the chirally symmetric phase even at zero temperature (even though the theory is still 
asymptotically free for A^f not too much larger than A^^ ''). The shape of the phase boundary 
near the critical flavor number, A^f < N^"^ , is most interesting from our viewpoint. In this 
region, the critical temperature is very small, and thus the system is probed in the deep 
IR. As a main result of this paper, we have shown that this connection becomes most 
obvious in an intriguing relation between the shape of the phase boundary for A^f < A^" 
and the IR critical exponent 6 of the running coupling at zero temperature. In particular, 
the flattening of the phase boundary in this regime is a direct consequence of |6| being 
smaller than 1. Since both the shape of the phase boundary and the critical exponent are 
universal quantities, their relation is a generic prediction of our analysis. It can directly be 
tested by other nonperturbative methods, even though it may numerically be expensive, 
e.g., in lattice simulations. 

Let us now critically assess the reliability of our results. Truncating the effective action, 
at first sight, is an uncontrolled approximation which can a priori be justified only with 
some insight into the physical mechanisms. The truncation in the quark sector support- 
ing potential critical dynamics is an obvious example for this. The approximation can 
become (more) controlled if the inclusion of higher-order operators does not lead to seri- 
ous modifications of the results. In the quark sector, it can indeed easily be verified that 
the contribution of many higher-order operators such as (^^)^ or mixed gluonic-fermionic 
operators is generically suppressed by the one-loop structure of the flow equation or the 
fixed-point argument given below Eq. (|T7jl . This holds at least in the symmetric regime, 
which is sufficient to trace out the phase boundary. By contrast, we are not aware of similar 
arguments for the gluonic sector; here, higher-order expansions involving, e.g., {F^yF^^Y 
or operators with covariant derivatives or ghost fields eventually have to be used to verify 
the expansion scheme. At finite temperature, the difference between so-called electric and 
magnetic sectors can become important, as mediated by operators involving the heat-bath 
four- velocity m^, e.g., {F^^UuY . view of results obtained in the Landau gauge [SHI, the 
inclusion of ghost contributions in the gauge sector appears important if not mandatory for 
a description of color confinement. A posteriori, the truncation can be verified by a direct 
comparison with lattice results. In the present case, this cross-check shows satisfactory 
agreement. 

The stability of the present results can also be studied by varying the regulator. Since 
universal quantities are independent of the regulator in the exact theory, any such regulator 
dependence of the truncated system is a measure for the reliability of the truncation. As 
was already quantitatively verified at vanishing temperature in |5n| . the present quark 
sector shows surprisingly little dependence on the regulator which strongly supports the 
truncation. By contrast, we do not expect such a regulator independence to hold in the 
truncated gluonic sector. If so, it is advisable to improve results for universal quantities 
towards their physical values. This can indeed be done by using stability criteria for the 
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flow equation which have lead to optimization schemes [HHl O EHl ■ We expect that the use 
of such optimized regulators give better results for dimensionless quantities, e.g. Eq. (|i9|l 
or the IR critical exponent 9. In any case, we have confirmed that, for instance, the linear 
regulator [SHI, which satisfies optimization criteria in various systems, leads to the same 
qualitative results as presented above. Further regulator studies are left to future work. 

Further generalizations of our work will aim at a quantitative study of the effect of finite 
quark masses; the formalism of which has largely been developed already in this work. 
Owing to the mechanism of fermionic decoupling, we expect that the largest modifications 
arise from a realistic strange quark mass which is of the order of the characteristic scales 
such as or the scale of xSB. 

Let us finally stress that our whole quantitative analysis relies on only one physical 
input parameter, namely the value of the gauge coupling at a physical input scale. This 
clearly demonstrates the predictive power of the functional RG approach for full QCD, and 
serves as a promising starting point for further phenomenological applications. 
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A Thermal moments and threshold functions 
A.l Thermal moments 

Let us first define the auxiliary functions / which are first introduced in Eq. ((221): 




(A.l) 




(A.2) 



r(«) 



1 1 



■u coth(M) 



(A.3) 




(A.4) 






1 1 u 

2 u'^d sinh u ' 



(A.5) 




2V^v V / dx{27ivxY^^T{-ed,{27ivxYu)cos{27rqx). (A.7) 

g=-oo 
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Here, the sum over q arises from the apphcation of Poisson's Formula to the (usual) Mat- 
subara sum. These functions are needed for the construction of the thermal moments Jij 
hf, gf, Hf and Gf which are related to the regulator via Eqs. © and ((121) , by 



:= h'^[m,v) = 

J 


poo 

/ ds h"^ {3,771)3^ ff(^s,v) , 
'0 


(A.8) 


hf:=hf{v) = 


poo ^ 

\ dsh{s)s^f^{s,v) , 
Jo 


(A.9) 


9f--=9t{v) = 


POO 

/ dsg{s)s^fT{s,v), 
Jo 


(A.IO) 


Hf:=Hf{v) = 


POO ^ 

/ dshis)s'f^{s,v), 
Jo 


(A.ll) 


Gf:=Gf{v) = 


POO 

/ dsg{s)s^ ff{s,v) , 
Jo 


(A.12) 



where m denotes a dimensionless quark mass parameter. It is more convenient to express 
the moments in terms of the regulator functions h{y) and g{y) in momentum space which 
are defined in Eqs. © and (0). In order to obtain the representations for hj and gj, we 
introduce 



s 



b+l roo 



duuV'" = l (6>-l) (A. 13) 



r(fe + l) 

in Eq. (IA.9I1 . (lA.inil and (IA.8I1 and use Eq. (lA.lll and (IA.2I1 . respectively: 



q=—oo 



y=0 

(A.14) 



(A.15) 



= V{h + l)V^ }-J, ^--H^^j \-Ty) L duug{y + u + x) 

Note that b is an arbritrary parameter which can, e.g., be used to avoid fractional deriva- 
tives. By applying Poisson's formula to the (usual) Matsubara sum, we have obtained the 
sum over q which converges fast for k >T. Moreover, we need Hf and Gf for j = 0, which 



(A.16) 
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are used in App. [B1 Integrating Eq. ()A.11|1 and (jA.12|l by parts and using Eq. (fin and 
()A.7j) . we obtain 



Hn. = —— — — > ax cos 0— / du Kii[y + u + x ) 



y=0 

(A.17) 



2 / X\ ( d\b-^^ . 

Gn = —r, ; — 1= 7 dx COS 0— / du -oiy + U + X 

° r{b + l)^ ^£^^Jo Vv) \ dyJ io u + x^'^' 

(A.18) 

In this paper, we use the exponential regulator. For the gluon and ghost fields, this 
regulator is given by 

i?,(A) = Ar(A) with r(y) = -^, (A.19) 

gy — 1 

and the functions h{y) and g{y) read ^ 



^^y^ = ^iz^ 9(y) = ^~'- (A-20) 



For the quark fields, the exponential regulator reads 



Rtm = ^Pr^{^) with rM = ^=L=-l, (A.21) 



a/1 - e~y 

and the functions h'^i^y, ^) and g^{y, ^) are given by 



h^(y,m) = — and q^(y,m) = — ^ . 

(A.22) 

Inserting Eq. (IA.2fljl and (1A.22I1 into Eqs. ((A.14jl - (1A.18)1 completely determines the desired 
thermal moments. 

A. 2 Threshold functions 

In Sec. m the regulator dependence of the fiow equations of the four-fermion interactions 
is controlled by threshold functions. The purely fermionic threshold functions are defined 
by 

where t = T/k and w are dimensionless quantities, the latter being associated with finite 
quark masses. Dots denote derivatives with respect to y^. The dimensionless momentum 
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= + y depends on the (dimensionless) fermionic Matsubara frequencies Un = (2n + 
l)7rt. The function p^{y^) is related to the regulator shape function by 



r^{y^) = ^l^^^^-l. (A.24) 
The factor v^^ is proportional to the volume of the d dimensional unit ball: 



= 2'^+Vir (0 . (A.25) 

In Sec. m we only need Using the exponential regulator Eq. (jA.2Hl and u; = for 
massless quarks, the fermionic threshold function l[^\t,0) reads, 

oo 

lf\t,0)= ^(-ire-t, /f)(t^o,0)^l. (A.26) 

n=— oo 

The threshold functions l^^l'''{t, Wi, W2) arise from Feynman graphs, incorporating fermionic 
and bosonic fields: 



„, ^ poo 

Cnr(^,^i,^2) = -^tJ2 dyy 



[p^{y-ii?i + wiY^[pA{yA) + 



f rii[p^{y^) - y^p4y^)] ^ n2[pA{yA) - yAPAivA)] \ .^27) 
I pAv^) + ^1 Pa^Va) + W2 J ■ 

Here, wi and W2 are dimensionless arguments and dots denote derivatives with respect 
to y^ and y^, respectively. In analogy to the fermionic case, the dimensionless bosonic 
momentum = + 2/ depends on the (dimensionless) bosonic Matsubara frequencies 
= An'^n^t^ . The (bosonic) regulator shape function r is connected with pa by the 
relation 

PA{yA)=yA[l+r{yA)]. (A.28) 

In Sect, m we need /[^^^^ and l^i-f^^. Using the exponential regulator Eq. (lA.lQjl and (IA.2H1 . 
we can calculate the integrals analytically in the limit t — ^ and Wi = W2 = ioY d = 4: 

\imlff^\t, 0, 0) = 1 and \iml%^^\t, 0, 0) = 31n(|) . (A.29) 

For t ^ 00 ov w ^ 00, the threshold functions l^^'^ and l^^^^ approach zero. For finite t 
and w, the threshold functions can easily be evaluated numerically. 
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B Resummation of the anomalous dimension 



Here, we present details for the resummation of the series expansion of the anomalous 
dimension rj, 



i=0 



The leading growth (l.g.) coefficients aj;f- read 

^iT{zd + m)r(m + 1) 



N{ 



22m _ 2 



r(^, + l) L ^ (2m)! 



r(2m) 2m-ed U ^ ^ (2m) ! ^ 2m-ed I fc ' fc ^'J ' l^-^J 

where i?2m are the Bernoulli numbers and is defined as 

z,:={d-l){Nl-l)c2. (B.3) 
The temperature and regulator-dependent functions Ci and C2 are given by 

c, = 2{H^{^)-Gt{^)), (B.4) 

Cl 

Note that Ci > and C2 > for |^ > 0. In the limits f ^ and |^ ^ oo, Ci and C2 are 
given by 

fc 

limc2 = c° = ^, (B.7) 

^lhn|ci = cr = 2v^(c(l + erf)-^), (B.8) 

fc 



T 



edC(l + ed) - 1 
2(C(l + e.)-^) 



hm C2 = = 2 , > (b-9) 



where we have used Eq. ()A.9j) - (|A.18j) for the exponential regulator and Q{x) denotes the 
Riemann Zeta function. 

Now, we perform the resummation of 77 along the lines of [iSI: We split the anomalous 
dimension Eq. into three contributions, 

V = Tlt + V2+^''. (B.IO) 
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where 77^ corresponds to the resummation of the term ~ T^B2m in Eq. ()B.2jl and 77^ to 
the resummation of the term containing the Nielsen-Olesen unstable mode (~ l/r(2m)) , 
representing the leading and subleading growth, respectively. The remaining contributions 
are contained in t]'^. 

First, we confine ourselves to SU{Nc = 2) for which the group theoretical factors are 
= iVc and = iV^ (1/4)™ = 2 (1/4)™ (see Appendix O for details), but we 
artificially retain the dependence in all terms in order to simplify the generalization to 
gauge groups of higher rank. 

We start with the resummation of rj^: For this purpose, we use the standard integral 
representation of the F functions [M] . 



T{zd + m)F(m + 1) 



dSi / dS2SiS2'^(.SiS2 



\m-l -{S1+S2) 



m—l 



where we have introduced the modified Bessel function 

Furthermore, we use the series representation of the Bernoulli numbers [HH] . 

B2m „(-!)""' ^ 1 



E 



\2m / ^ ^2m 
/=1 



(2m)! (27r)^ 

With the aid of Eq. (jB.llI) and (jB.13|) , we rewrite r/^ as follows 



(B.ll) 
(B.12) 

(B.13) 



7r2F(z, + l) 



m=l 1=1 



j 2Gpci 
\ 7r2/2 



V27r2/2 



m— 1-1 



In order to perform the summation over m, we define 

00 00 



(B.14) 



/2 ^ V/2 

=1 m=l 

00 00 00 



l=\ 771=0 n=— 00 

00 



fnx 



dxcosy — ) / c/tLii(e v^)af(x , t) 



(B.15) 



F(6+l)7g¥ _ ^ 

where we have used Eqs. (|A.9j) . (jA.lj) and (|A.13j) . The auxiliary function is defined as 



d \fc+3-ed 

dy) 



duu 



h{y+u+x^-t) + h{y+u+x^+t) 



(B.16) 



y=0 
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Using Eq. ()B.15jl . we obtain the final expression for 77^, 



^ A{d-2)N^G r, , . 



T(^,+l) 



which can straightforwardly be evaluated numerically. 

Now we turn to the calculation ofr]2, the subleading-growth part of r/. Here, a careful 
treatment of the zeroth Matsubara frequency which contains the Nielsen-Olesen mode, is 
necessary. More specifically, we transform the modified moments hf in Eq. (|A.15|) into 
a sum over Matsubara frequencies and insert a regulator function 7^(1-) for the unstable 
mode, 



4:71V / dsh{s)s-^e 



J„-sPn{v) 



(B.18) 



Here, we have introduced 



VJv) 



{2TTnv)'^ {n ^ 0) 
V(v) {n = 0) 



(B.19) 



The function V(v) specifies the regularization of the Nielsen-Olesen mode and is defined 
in Eq. (|D.3j) : the other modes with n 7^ remain unmodified. 
We rewrite 77^ by means of Eq. (jB.lHl , 



16N^G ^ 1 

r{z,+i) ^,r{2^) Jo 



A,reg /T^ 



2GpC] 



m—1 



(B.20) 



Now it is convenient to introduce an auxiliary function T^{q) which is defined as 



m=l 



rf2m) 



m— 1 _ 



g h 



A,reg 
2m-erf 



n^ + ^Kt-Jo Jo Jo 

, — 00 



(B.21) 



Here, we have used Eqs. ()B.18jl and ()A.13j) . Furthermore, we have defined the function 



— j / dt I duu^ 



^0 



h{y+u+ x—t^/^) + h{y+u+x +tv/^) 



(B.22) 
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Applying Eq. (|R2T]1 to Eq. dOHll . we obtain 



dpk^^^,{p)T,\2Gpc 



I5 fc J ) 



(B.23) 



which can straightforwardly be evaluated numerically. Finally, we have to calculate the 
contribution of the quarks to the gluon anomalous dimension. Performing analogous steps 
along the lines of the calculation of ?7j^, we obtain 



The auxiliary function Sf{q,m) is defined as 
1 



r(6+l)v^ 



E 



-1 



dx cos 



nx 



dtLiiie ) erf (m, x^, t, m) 



where af{u,x,t,m) is given by 



(B.25) 



crf{u, X, t, m) 



d \b+3-ea 

dy) 



duu 



hf (y/y+u+x-t, m] + hU^/y+u+x+i 



, m 



+ht ( -^yy+u+x-t,m) + h'^A-\fy+u+x+t,m 



j/=0 



(B.26) 



The regulator function occurs in the function hf{^yy,m) which is related to h'^{y,m) by 

htiVy^m)^hAy,m). (B.27) 
There is one essential difference between the resummation of 77^2 that of v'^: the 



regulator shape function r{y) can be expanded in powers of y, while the corresponding 
function r^{y) for the quark fields should have a power series in ^ which is a consequence 
of chiral symmetry [HII; this explains the notation hf{^yy,fh). 

We stress that all integral representations in Eqs. (|B.17j) . ()B.23j) and ()B.24jl are finite 
and can be evaluated numerically. For d = 4 and in the limit T ^ 0, the results agree 
with those of Ref. [13]. 

The remainder of this section deals with a generalization to higher gauge groups. Since 
we do not have the explicit representation of the color factors Tm for gauge groups with 
A'c > 3 at hand, we have to scan the Cartan subalgebra for the extremal values of and 
r^. However, as discussed in App.lOl these extremal values of and can be calculated 
straightforwardly. Their insertion into Eq. (jB.2j) allows to display the anomalous dimension 
for SU{3) in terms of the already calculated formulas for SU{2): 



„SU{3) 

'/s 



„SU(3) 

'/s 



A^c^S 3 



Afc-^3,ci-»3ci/4 9 



H — TJ 

Nc^3,ci^ci/A 3 

4 



H — rj 

7Vc-»3,ci^ci/3 9 



Afc-^3,ci-»4ci/3 



(B.28) 
. (B.29) 
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The notation here serves as a recipe for replacing N^. and ci, defined in Eq. ()B.4jl . which 
appear on the right-hand sides of Eqs. (fRTTll . (jR23ll and (jR24jl . Note that the replace- 
ment of Nf. results also in a modification of z^, defined in Eq. ()B.3|1 . However, C2, which 
appears in the definition of Zd, remains unchanged for all gauge groups and depends only 
on the dimension d. 

C Color factors 

In the following, we discuss the color factors and rf which carry the information of the 
underlying SU{Nc) gauge group. First, we summarize the discussion of Ref. [IHlIiniEni for 
the "gluonic" factors r/^ appearing in the flow equation: Gauge group information enters 
the flow of the coupling via color traces over products of field strength tensors and gauge 
potentials. For our calculation, it suffices to consider a pseudo-abelian background field A 
which points into a constant color direction Therefore, the color traces reduce to 

n'^in"^ . . . tv^lT^^'T"^ . . . T"^')] , (C.l) 

where the parentheses at the color indices denote symmetrization. These factors are not 
independent of the direction of n"-, but the left-hand side of the flow equation is, since it 
is a function of the n"-independent quantity \F^j^F^^. For this reason, we only need that 
part of the symmetric invariant tensor trc[T*^"^ . . .T'^^^)] which is proportional to the trivial 
one, 

trjT^'^^T^^ . . . T'^-)] = r, . . . + . . . . (C.2) 

Here, we have neglected further nontrivial symmetric invariant tensors, since they do not 
contribute to the flow of Wfc(^), but to that of other operators which do not belong to our 
truncation. For the gauge group SU(2), there are no further symmetric invariant tensors 
in Eq. ()C.2jl . implying 

^su(2)^2, z = l,2,... . (C.3) 

However, for higher gauge groups, the above mentioned complications arise. Therefore, we 
do not evaluate the r/^'s from Eq. (|C.2)1 directly; instead, we use the fact that the color 
unit vector can always be rotated into the Cartan sub-algebra. Here, we choose the 
two color vectors n"' which give the extremal values for the whole trace of Eq. (jC'.ljl . For 
SU(3), these extremal choices are given by vectors ra" pointing into the 3- and 8-direction 
in color space, respectively: 

r5'""''=2+j^. -r"'"=3gy". (C.4) 

Finally, we turn to the color factors of the quark sector. The above considerations also 
hold for the contributions of the flow equation which arise from the fermionic part of our 
truncation Eq. ([T5|) and (fT7|l . Taking into account that quarks live in the fundamental 
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Figure 8: Dependence of the critical temperature on the smeared regularization of the 
Nielsen-Olesen mode with m labehng the regulator. The left and the right panel show the 
results for N^. with A^f = 3 and A^f = 11 massless quark flavors, respectively. The limit 
m — i> oo can be identified with the stationary point, and thus optimal regulator, in the 
class of considered regulators. This justifies constructively the procedure used in the main 
text which was derived from general considerations. 



representation and choosing a color vector pointing into the 3- or 8-direction, we obtain 

r*r'" = 2[^\ '■5'"'"= 2(1^) '+(5)* ^ = 1.2..... (C.5) 

Again, all complications are absent for SU(2) and we find r/"'^^^^^-* = rff^^^K 

The uncertainty introduced by the artificial n°- dependence of the color factors is the 
reason for the uncertainties of our results for the critical temperature and the fixed point 
values in three and four dimensions. 



D Regulator dependence from the unstable mode 

In this section, we discuss the regulator dependence of the critical temperature T^, arising 
from the details of projecting out the unstable Nielsen-Olesen mode. As already explained 
in the main text, removing the tachyonic part of the unstable mode corresponds to an exact 
operation on the space of admissible stable background fields. In the present context, it 
even suffices to remove only the thermal excitations of the tachyonic part of the mode, since 
the imaginary part arising from quantum fiuctuations can easily be identified and dropped. 
In the following, we take a less strict viewpoint and allow for a smeared regularization of 
this mode in a whole class of regulators. 

Since the true physical result will not depend on this part of the regularization, we can 
identify the optimal (truncated) result with a stationary point in the space of regulators, 
using the "principle of minimum sensitivity", cf. In order to inhibit the thermal 

population of the Nielsen-Olesen mode E^*^ at finite temperature, it suffices to regularize 
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only the soft part (zero Matsubara frequency) of this mode as follows: 



fc^ ^ k 

The function V{^) has to satisfy the following constraints 



va) + ^^\^. (D.i) 



limP(^) = and lim P(^)^oo. (D.2) 

T/fc->0 " T/fc-oo " 



In the following, we choose 

P(f ) = Vmij) = ijr with m > (D.3) 

as a convenient example. As a regulator optimization condition, we demand that T^^ should 
be stationary with respect to a variation of the optimal regulator function. Calculating 
as a function of the parameter m, the optimization condition for the regulator function 
translates into 



dT,r 



dm 



0. (D.4) 



The solution m = m defines the desired optimized regulator. 

As an example, we show Tcr{m) /Tcr (oo) as a function of m for = 3 with iVf = 3 
and with A^f = 11 quark flavors in Fig. |H1 We find that the optimized regulator is given by 
m ^ oo for all N^. and Af. This represents an independent and constructive justification 
of the regularization used in the main text, corresponding to the choice m oo. 
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